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We report a series of Lattice-Boltzmann simulations of the sedimentation velocity of charged 
, disks. In these simulations, we explicitly account for the hydrodynamic and electrostatic forces on 

■ disks and on their electrical double layer. 

By comparing our results with those for spheres with equal surface and charge, we can clarify the 
(N ■ effect of the particle shape on the sedimentation process. We find that disks and spheres exhibit a 

different dependence of the sedimentation velocity on the Debye screening length. An analysis of the 
behavior of highly charged disks (beyond the scope of the linearized Poisson-Boltzmann equation) 
shows that, in that regime, the charge dependence of the sedimentation velocity of disks and spheres 
^-j- . is similar. This suggests that, at high charge, the effective hydrodynamic shape of the disks becomes 

' more spherical. 

I 1 PACS numbers: 82.70.Dd, 82.70.-y, 66.10.-x, 47.65.+a, 82.45.-h 

o 

& • I. INTRODUCTION 

a 

' Suspensions of charged disks are of great practical importance. Examples range from clay suspensions to blood. In 
the present paper, we present calculations of the electrokinetic behavior of charged disks. As disks are not spherically 
symmetric, they also provide an ideal model system to study the effect of shape on the coupling between electrostatic 
O , and hydrodynamic response of a macroscopic particle. 

There exists an extensive experimental literature on the transport properties of suspensions of charged disks. Yet, in 
y—{ \ spite of the importance of these systems, there is surprisingly little theoretical knowledge about the effect of the charge 
. of the disks on their transport properties. One reason may be that the non-spherical geometry greatly complicates 
the use of the analytical approach that is used to describe charged, spherical particles. Whilst there are papers that 
consider the hydrodynamical properties of uncharged disks or the electrostatic properties of charged disks 0, 0] , 
we are not aware of any theoretical publications that treat the interplay between electrostatics and hydrodynamics 
for charged disks. 

Several theoretical studies suggest that, in general, there may be non-trivial coupling effects due to shape asymme- 
tries 0, 0| . The need for a numerical (rather than an analytical) approach is related to the fact that the analytical 
approaches are usually only tractable in certain limits. For example, to arrive at tractable analytical expressions it is 
often necessary to assume that the Poisson-Boltzmann equation can be linearized or that the Debye screening length 
is small compared to the linear dimensions of the charged colloid. Yet these conditions are rarely satisfied for charged 
biocolloids and biopolymers. These (non-spherical) particles usually have high charges under physiological conditions 
making the use of the linearized Poisson-Boltzmann equation questionable. 

The total number of ionizable groups of colloids can be determined experimentally by titration. However, the 
experimental determination of the effective charge of colloids usually exploits non-equilibrium techniques, such as 
• • ■ electrophoresis or sedimentation. The interpretation of many of these experiments is based on theoretical expressions 
. ' that have been derived for weakly charged, spherical colloids. Yet, it is not at all obvious that expressions valid for 
spheres can be extrapolated to disks. Hence, it is important to understand the effect of shape and charge on the 
$H ' transport behavior of non-spherical charged particles. 

In the present paper we report a numerical study of the electrokinetics of charged disks. In order to perform these 
simulations, we extended a Lattice-Boltzmann scheme that we had developed to study the electrokinetic behavior 
of spherical colloids. With this simulation method, we can study the dynamics of charged particles at the Poisson- 
Boltzmann level. Hence our simulations ignore charge fluctuations and static charge correlation effects. In the absence 
of polyvalent ions, the Poisson-Boltzmann description is expected to work quite well. 
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Our simulations extend from the regime where the electrohydrodynamic coupling is small (low charge, thin double 
layer) to the case where the charge is large and the double layer is extended. By performing simulations of many 
different conditions, we disentangle the effects of charge, double layer width and shape. 

The paper is organized as follows. In Sec.[n]we present the basic electrokinetic equations and the simulation method 
that we use to solve them. In order to study the effect of charge on the sedimentation of disks, we must first know 
the sedimentation behavior of uncharged disks. This is done in Sec. II I II where we use our model to compute the 
friction coefficients of a neutral sedimenting disk for motions parallel to and perpendicular to its axis of symmetry 
(we refer to these motions as "longitudinal" and "transverse", respectively). In Sec. IIVI we discuss the effect of the 
charge of the disk on the sedimentation velocity. In Sec. we consider the effect of the concentration of the disks 
on the sedimentation velocity. Sec. I VII focuses on the variation of the sedimentation velocity with Debye-screening 
length. Sec. IVIll contains a comparison of the sedimentation behavior of disks and spheres. Concluding remarks are 
contained in Sec. IVIIII 



II. ELECTROKINETIC MODEL 



We analyze a simple geometry in which one disk-like colloidal particle with radius a and height h sediments due to 
the action of a uniform external field. The disk has an overall charge Q = Ze, where Z is the valency and e is the 
elementary charge unit. The aspect ratio p is defined as p = 2a/ h . The disk is suspended in a symmetric electrolyte 
and, for the sake of simplicity, we assume that coions and counterions have the same mobility. The fluid mixture is 
characterized locally by the solvent density, p s , and by the microion electrolyte densities, p±. The latter also determine 
the local charge density of the fluid, q(r) = ze[p+(r) — p_(r)]. We restrict ourselves to monovalent electrolytes, i.e. 
z = l. 

On a macroscopic length scale, the dynamics of the system is governed by the standard electrokinetic equations [f| 
that specify the interplay between the electrical potential, local charge density, electrical currents and fluid flow: 



dt 



Pk 



d k T 
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(1) 
(2) 



jfc = -Pfcv + D k [Vp k + z k p k V$] 



(3) 



where rj is the shear viscosity, P is the pressure, v is the fluid velocity, fcgT = 1 measures the temperature, and 
D k stands for the diffusivity of each electrolyte species (which reduce to a single constant for symmetric electrolytes). 
$ is the electrostatic potential, while <& = $(/csT/e) is an appropriate dimensionless potential which satisfies the 
Poisson equation 



k=± 



(4) 



where lg — (3e 2 / (Aire) is the Bjerrum length, e = t^e r denotes the dielectric constant of the medium, whilst p w refers 
to the charge density due to embedded solid objects, either colloids of solid walls. 

Equation Q simply expresses a conservation law and Eqn. is the constitutive equation. Together with the 
incompressibility condition V • v = 0, Eqn. corresponds to the Navier-Stokes equations for an incompressible, 
isothermal electrolyte. In the presence of external forces (such as the gravitational field), the corresponding force 
must be added to the right-hand side of Eqn. 



A. Simulation method 



To simulate the sedimentation of charged disks, we used the lattice-Boltzmann scheme reported in Ref. |2|. We 
showed therein that lattice-Boltzmann can be used to compute transport coefficients of charged spherical colloids. 
Below, we briefly summarize the main features of the method and refer the reader to that reference for further details. 

The Lattice-Boltzmann (LB) method is the lattice counterpart of the Boltzmann equation. It prescribes a dynamical 
evolution rule for the distribution function rij(r, i), which represents the density of particles at the lattice node r, at 
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the discrete time t and with the discrete velocity Cj . The density- weighted moments of the local velocity distribution 
correspond to the hydrodynamic fields. In particular, J^Cjn^r, t) = pv is the fluid's momentum; it satisfies the 
Navier-Stokes equation on length and time scales that are large compared to the lattice spacing and the LB time step, 
respectively. 

The electrolyte species are simulated by following the diffusion and convection of the local densities of coions 
and counterions described by Eqns. and JSJ- This equation is based on the flux of each species along the links 
that connect neighboring nodes, and ensures strict local charge conservation. This local charge, combined with the 
corresponding electrostatic potential — computed by a numerical solution of the Poisson equation J1J — provides the 
local force that accelerates the fluid. With this technique, colloidal particles are simply introduced as surfaces where 
the collision rules of the populations of the neighboring nodes are modified to ensure non-slip boundary conditions ||| . 
The link-based definition of the flux of the electrolyte species leads to a straightforward implementation of the no-flux 
boundary condition for each of the ionic species at solid surfaces. This suppresses possible charge leakage through the 
solid walls. 

For reasons of computational convenience, we choose the value of the kinematic viscosity v = 1/6 (in lattice units) @ 
and the p s = 1, as the density unity. The external (gravitional) field that induces sedimentation was chosen to be 
1CP 6 , a value that is well inside the linear-response regime. This gravitational field generates fluid velocities of the 
order of 1CP 8 (again, all expressed in lattice units). The diffusivity of the electrolyte is set to D = 0.19, a value for 
which spurious diffusion due to lattice advection (see Ref. 0) is negligible. The values of the diffusivity and the flow 
velocity correspond to Peclet numbers smaller than 1CP 1 . In the simulations described in the subsequent sections we 
vary the salt concentration between 7 x 1CP 4 and 5 x 1CP 3 as a way to control the electrical double layer thickness. 

III. SEDIMENTATION OF NEUTRAL DISKS 

Before assessing the role of electrostatics on the sedimentation of non-spherical particles, we performed LB simula- 
tions to compute the sedimentation velocity of uncharged hard disks. Such reference calculations are needed because, 
in contrast to the case of hard spheres, analytic expressions for the sedimentation velocity of an isolated hard disk 
only exist in the limit of infinitely thin disks j^, E3 ■ We are not aware of analytical results for disks with finite aspect 
ratios. 

We have simulated the sedimentation of disks with two different nominal aspect ratios p = 10 and p = 5, corre- 
sponding to disks of lateral dimension h = 2, and radii a = 10 and 5, respectively. However, these aspect ratios 
are only approximate: in the LB approach, the hydrodynamic boundary of a solid particle is usually located close 
to the midpoint of links joining fluid and solid nodes. In practice, the hydrodynamic shape of an object may differ 
slightly from the nominal one. For this reason, we need to calibrate the shapes of the disks. An unambiguous way 
to determine these effective sizes (the hydrodynamic radius and height) would be to measure the friction coefficients 
of the particles and compare the results with the appropriate analytical expression of an infinitely thin disk. Great 
care must be taken when comparing the LB numerical results for the friction coefficient with results obtained for an 
isolated disk. Since we use periodic boundary conditions, the simulations measure the friction coefficient of a regular 
array of particles at volume fraction ip = ira 2 h/L 3 , where L is the diameter of the simulation box. Hence, in order to 
extrapolate to infinite dilution, we must perform a series of simulations with increasing box size. The same procedure 
will also be used to compute the sedimentation velocity of a charged disk in the dilute limit. 

A. Friction coefficients 

In the following calculations, the reference frame is centered on the disk and the disk is fixed on the lattice. 
Hence, the gravitational force acting on the disk becomes a body force, with opposite direction, that acts to the 
fluid. We then determined the total fluid velocity at steady state, vp. As this velocity is equal and opposite to the 
sedimentation velocity of the particle, Ud = — Vf, we can relate the sedimentation velocity to the friction coefficient 
through Ud = Fg/t;, where Ud is the sedimentation velocity of the disk, £ is its friction coefficient, and F g is the 
applied gravitational force. We computed both the friction coefficient for motion along the symmetry axis of the disk 
and the one for perpendicular (transverse) motion. 

Hashimoto [ll| has shown that the friction coefficient of an array of hard spheres depends, at low volume fractions 
ip, as (£/£o) -1 = 1 — 1 .76C/3 1 / 3 + <p + 0(p 2 ), which shows that the initial decrease in the sedimentation velocity of an 
array of hard spheres is controlled by the colloid volume fraction. 

In Figure ^ we show the transverse and longitudinal friction coefficients for a neutral disk with p — 5 at various 
volume fractions, normalized by the Stokes friction coefficient of a sphere with the same area, i.e. by £a = Gt^t/R, 
R = \Ja(a + h)/2. The figure shows that, just as in the case of neutral spheres, the friction varies linearly with p 1 / 3 . 
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FIG. 1: Transverse (circles) and longitudinal (squares) friction coefficients for a sedimenting neutral disk of aspect ratio p — 5 
as function of the volume fraction of the array. The friction coefficients are normalized by the friction coefficient of a sphere 
with the same radius as the cylinder, i.e. £a = (SurjR. The dashed lines are linear fits. 



Knowledge of this concentration dependence allows us to extrapolate the numerical results to estimate the friction 
coefficients and £m) at infinite dilution. One limiting case is known: the friction coefficients of an infinitely thin 
disks (p — > oo) is identical to that of an oblate spheroid with the same aspect ratios between its main axis |10|. 

We summarize the results for the normalized friction coefficients for disks of two aspect ratios in Table HJ For the 
disk with p = 5, we conclude that, to a good approximation, the hydrodynamic radius and height correspond to the 
nominal ones. The numerical values obtained for the larger disk (p — 10) are less satisfactory than for the shorter 
one. While one could conclude that, for a larger aspect ratio the hydrodynamic radius and height are different with 
respect to the nominal ones, and hence recompute the effective values for "a" and h so as to adjust £a, this does not 
seem satisfactory, since for a larger object the disagreement between the nominal and hydrodynamic size is expected 
to decrease. Moreover, because we have always used cubic simulation boxes, we can attribute such a difference to 
the stronger coupling between image disks for longitudinal sedimentation. On the other hand, the ratio between 
the perpendicular and parallel frictions do agree with the value estimated on the basis of the nominal size, 

suggesting that the deviations come mostly from uncertainties related to the extrapolation from finite volume fraction 
values. Hence, we conclude that, also for this shape, the disagreement between the two sizes is negligible, and we 
ascribe the deviations to interactions with the periodic images. 



IV. SEDIMENTATION VELOCITIES OF CHARGED DISKS: CHARGE DEPENDENCE 



Before discussing the computed sedimentation velocities of charged disks, we briefly recall the theoretical results 
concerning the sedimentation velocity of weakly charged spheres, since this theory serves as a reference point for the 
discussion of the results for disks. Booth ^| (and Ohshima et al. predicted that the sedimentation velocity 

U S (Z), where Z is the valence, of an isolated sphere is a quadratic function of the sphere valence, which can be 
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£a = 6irriR 


p=5 


p=10 




Simulation 


Oblate spheroid 


Simulation 


Oblate spheroid 




0.91 


0.9 


0.93 


0.85 




1.1 


1.06 


1.19 


1.11 


Cx/C|| 


0.83 


0.84 


0.78 


0.77 



TABLE I: Friction coefficients of an isolated disk in transverse and in longitudinal motion normalized by the Stokes friction 
coefficient £a = 6ivriR of a sphere with equal surface area, i.e. with radius R — */ a(a + h) /2 (a is the radius of the cylinder). 
We compare the computed value of the friction coefficients with approximate theoretical values for an oblate spheroid with the 
two axis equal to the disk radius and to the disk height. We studied two disks with aspect ratios p = 5 and p — 10, respectively. 



expressed as 

Us{Z) =1-c 2 (kR)Z\ (5) 



U s (0) 

where U s (0) corresponds to the sedimentation velocity of a hard sphere in the dilute limit. In the regime where the 
Debye-Hiickel theory is valid, the pre- factor C2 for a symmetric 1-1 electrolyte is given by 

= rt^D f ^ (6) 

where R is the radius of the sphere and, with p° denoting the bulk value of the charge density of species k, k = 
^/47tTb JTu z 'kPk stands for the inverse Debye length that characterizes the size of the electrical double layer. The 
function f(nR) is defined as 

1 

M) 2 



f( K R) = 1__ \ g 2kR (3E 4 (kR) - 5E 6 {nR)) 2 + 8e KR {E 3 {nR) - E 5 {nR)) 



(4:E 3 (2kR) + 3E 4 (2kR) - 7E 8 (2kR))] , (7) 



expressed as a linear combination of the integral function, E n (a;) = x 11 ^ 1 dt t~ n exp(— t). 

It is reasonable to assume that, in the case of disks, the dependence of the sedimentation velocity on colloidal charge 
in the Debye-Hiickel limit has the same functional form as Eqn. JSJl, where all the shape dependence and hydrodynamic 
coupling enters through the factor C2(kR). To test this, and to analyze the role of charge on the sedimentation velocity 
of disks, we performed a series of simulations at constant Debye screening length and volume fraction (ip = 7.2 x 10~ 4 
for p — 5, and (p = 2.9 x 1CP 3 for p — 10). Since in many cases of practical interest colloidal particles are highly 
charged (see e.g |l ll for the case of disk-like clay particles) we performed numerical simulations covering a wide range 
of disk charges. 

The results obtained are displayed in Fig. [21 where the velocity, normalized by the sedimentation velocity of a hard 
disk at the same volume fraction, is depicted as a function of the surface charge a = eZ/ [2na(a + h)]. One can 
identify the quadratic dependence at low charge, consistent with Booth theory for spheres. Such a dependence can 
be clearly appreciated in the figure's inset. After this quadratic growth, a crossover region is identified, for surfaces 
charge densities between 0.1 and 0.4, before entering the asymptotic regime at even larger surface charge densities, 
where the sedimentation velocity increases much more slowly. This behavior is consistent with numerical results 
on the sedimentation velocity of charged colloidal spheres which also show a deviation from Booth's predictions for 
surface charge densities around 0.1. 

Comparing the deviation of the sedimentation velocity for transverse and longitudinal motion, as displayed in 
Fig. [3 one can see that the sedimentation velocity decreases faster with charge for transverse motion, regardless of 
the aspect ratio. This is indicative of a stronger electrokinetic coupling for transverse motion. Interestingly, the 
sedimentation velocities (expressed by 1 — Ud(Z)/Ud(0)) become almost independent of the aspect ratio of the disks 
for large surface-charge densities. In the same figure, we also plot the decrease in sedimentation velocity for a sphere 
of radius R — 4.2, which shows the same dependence on surface charge as that of the disks. 

Figure [3 allows us to draw some qualitative conclusions concerning the nature of the errors that are made when 
estimating the charge of disk-like particles by assuming the validity of Booth's theory 01 . As mentioned above, 
Booth's theory is valid only for weakly charged sedimenting spheres. As we now have numerical results for the 
sedimentation velocity of disks, we can identify two sources of errors. At low charge, where the quadratic dependence 



6 




FIG. 2: Surface-charge-density dependence of l — Ud(Z)/Ud(0) for a disk with aspect ratio p — 5 (squares) and p = 10 (triangles) 
in transverse (dotted lines) and longitudinal (continuum lines) motion at ka = 1. For comparison, we also show the behavior 
for a sphere with radius R = 4.2 (circles). In the inset, we show 1 — Ud{Z)/Ud(0) for weakly charged disks as function of Z 2 . 
The curves for the disk with p = 10 are plotted against Z 2 /10 to show them on the same scale. The inset shows that for small 
charges, the sedimentation velocity varies quadratically with charge. The lines in the inset are linear fits to the simulation 
data. In contrast, the curves in the main figure are simply meant as a guide to the eye. 

of the sedimentation velocity on colloidal charge holds, there will still be some discrepancy in the coefficient c<x. We 
address this issue in more detail in Section lYlII At high charge, where even the quadratic surface-charge dependence 
does not apply, the use the Booth theory leads to serious errors in the estimation of the zeta potential. To give an 
idea of the magnitude of this error, one should compare the curves shown in Fig. |2 with a parabolic extrapolation 
of the curves up to a ~ 0.1. Our calculations also suggest that if the sedimentation velocity of a charged colloid is 
plotted as a function of the surface-charge density (and not as a function of the total charge as one might be tempted 
to do when comparing with Booth's theory) the curves corresponding to the sedimentation velocity of a variety of 
charged disks show the same functional dependence. 

V. SEDIMENTATION VELOCITIES OF CHARGED DISKS: VOLUME FRACTION DEPENDENCE 

As discussed above, the inverse friction coefficient of a dilute, ordered array of hard spheres scales as tp 1 / 3 . In 
Ref. pi, we have verified that this functional dependence also holds for charged spheres, provided that the system is 
dilute enough to guarantee that there is no significant overlap of the double layers. For charged disks, we expect the 
ip 1 / 3 dependence to hold under the same circumstances. In order to test whether there is a detectable effect of the 
overlap of electric double layers of different disks, we have computed the normalized friction coefficients for disks of 
aspect ratio p = 5 as a function of the volume fraction, for volume fractions up to 10%, for different widths of the 
diffuse layer. 

In Fig.[3]we show the results for a weakly charged disk, both for transverse and longitudinal motion. Note that, in 
the dilute limit, the friction coefficients will depend on Ka due to the electrohydrodynamic interaction. For transverse 
motion the convergence to the dilute limit is slower, indicating a stronger coupling between disks; we attribute this to 
the fact that the distance of closest approach coincides with the external field direction. Although for na — 1/2 and 
high volume fractions the diffuse layers overlap, the effect of the diffuse layer is much weaker than the volume-fraction 
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FIG. 3: Volume-fraction-dependent normalized friction coefficients for a disk with radius a = 5, aspect ratio p — 5, and valency 
Z = 10 for na equal to 0.5 (spheres), 0.8 (pluses), and 2.1(filled diamonds). The upper curves are for the transverse friction 
coefficient, while the lower curves are for a longitudinal friction coefficient. The curves are normalized by the friction coefficient 
of an isolated sphere with equal surface area = &nr)yja(a + h)/2 (see Table for the correspondent neutral values). Curves 
are drawn as a guide to the eye. 



dependence or than the effect of the particle shape. The shape effect is reflected in the substantial difference between 
the friction coefficients for transverse and longitudinal sedimentation. Although barely visible in Figure there is a 
small but significant dependence of the friction coefficients on n. This we discuss in more detail in the next section. 

The dependence on k is better visible in Fig. where we display the friction coefficients of highly charged disks 
Although here the dependence on k is clearly visible, it becomes less important at higher volume fractions. This 
suggests that with increasing volume fraction, the effect of the overlap of diffuse layers becomes less important than 
the direct effect of hydrodynamic interaction between disks. At small volume fractions, we always recover the ip 1 / 3 - 
dependence of the inverse friction coefficients. However, deviations form this scaling are already noticeable at small ip. 
This fact indicates that the dilute regime is confined to smaller ip values for highly charged disks. In the limit na — > 00 
(vanishing Debye length), the sedimentation friction coefficient of a charged and a neutral disk should coincide. The 
fact that, for finite Ka, the ratio of the friction coefficients extrapolates at tp = to a value different from one indicates 
the importance of electrokinetic coupling for the sedimentation of charged disks. 

We stress that the values of k are computed on the basis of the electrolyte densities in the bulk. For concentrated 
suspensions of highly charged disks, the density of counterions added to the system to ensure charge neutrality 
(p_ = — Z/Vf, where Z is the valency of the sphere and Vf the volume occupied by the electrolyte) may exceed the 
concentration of added salt. In that case, the screening length in the salt "reservoir" is not simply related to 
the apparent screening length in the dense suspension. For our simulations this phenomenon becomes important only 
for the highest volume fractions (typically ip > 0.08). Hence, our extrapolation at low tp is not affected by this 
complication. 



VI. SEDIMENTATION VELOCITY OF CHARGED DISKS: DEBYE LENGTH DEPENDENCE 



Having analyzed the role of charge and volume fraction on the sedimentation velocity, we now consider in more detail 
the effect of the double layer width on the sedimentation of the disks. We follow the same procedure as in Section llll Al 
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FIG. 4: Volume-fraction-dependent normalized friction coefficients for a disk with radius a — 5 and aspect ratio p = 5 and 
valency Z = 100 or na equal to 0.5 (spheres), 0.8 (pluses), and 2.1(filled diamonds), a) Transverse friction coefficient; b) 
Longitudinal friction coefficient. In all cases,the curves are normalized by the friction coefficient of a sphere with equal surface 
area £a = Qi^rj^Ja{a + h)/2 (see Table U] for the correspondent neutral values). Curves are drawn as a guide to the eye. 



and study the sedimentation velocity Ud{Z) at different volume fractions, normalized by the corresponding velocity of 
isolated charged-neutral disks Ud(0). In this way, the ratio Ud(Z)/Ud(0) measures the reduction is the sedimentation 
velocity of one charged disk due to its electrokinetic interaction with the electrolyte. The reduction in sedimentation 
velocity in the dilute limit is interesting theoretically, because we can compare with analytic results for weakly 
charged spheres, although in experiments the reduced sedimentation velocity Ud(Z,ip)/Ud(0,tp) at finite ip is the 
relevant quantity. For simplicity, in the remaining part of this section, we will be writing Ud(Z) instead of Ud(Z,tp) 
but, unless explicitly stated, the volume fraction dependence is always assumed. 

In Figure EJa), we show the normalized sedimentation velocity as a function of the double layer width for a disk 
with aspect ratio p — 5 in transverse motion, for different volume fractions. For infinitely thin and infinitely broad 
diffuse layers, the sedimentation velocity should coincide with that of a charged-neutral disk, and hence the curve 
should approach one for both small and large na, as is indeed observed. The decrease at intermediate values of na is 
the result of the interplay between hydrodynamic dissipation and electrolyte diffusion. The largest effect is observed 
when the size of the double layer is of the order of the largest dimension of the disk, i.e. kcl ~ 1. The effect increases 
with decreasing volume fraction, consistent with the discussion in the previous section, and above volume fractions 
around 1% , the changes in normalized sedimentation become negligible. The minimum velocity also depends on 
volume fraction, an effect which is consistent with previous findings for spheres (la ]. A similar behavior is observed 
in Fig. [S^b), where longitudinal sedimentation for a weakly charged disk is depicted. It is interesting to note that the 
decrease in sedimentation velocity is slightly smaller. We can ascribe this effect to the fact that the distorted double 
layer is not isotropic and has a smaller contribution to the friction when the wider side of the disk is exposed to a 
region where the velocity gradients are smaller. 

In Figures EJa) andHJb), we show the sedimentation velocity for a weakly charged disk with a somewhat larger 
aspect ratio: p = 10. The trends are the same as for the disk with p = 5, although the minimum velocity seems to 
depend on aspect ratio, and is achieved now for slightly narrower double layers. The reduction in absolute terms is 
now smaller, but this is simply due to the lower surface charge density as compared with the smaller {p = 5) disk. It 
is worth mentioning that for disks with p = 5 we could effectively reach the dilute limit for Ud(Z)/Ud(0). In contrast, 
for disks with p — 10 we had to perform the dilute-limit extrapolation. 



9 




(a) Transverse motion 




(b) Longitudinal motion 



FIG. 5: Normalized sedimentation velocity of a charged disk with aspect ratio p = 5, radius a = 5, and valency Z — 10 as 
function of the EDL thickness expressed in dimensionless units so. The different simulation points correspond to ip = 3.8 x 1CP 2 
(circles), 4.8 x 1CT 3 (squares), 7.3 x 1CP 4 (pluses), 4.6 x 1CP 4 (stars), and 3.1 x 1CF 4 (diamonds), (a) Transverse motion; (b) 
Longitudinal motion. Curves are drawn as a guide to the eye. 
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FIG. 6: Normalized sedimentation velocity of a charged disk with aspect ratio p = 10, radius a = 10, and valency Z = 10 as 
function of the EDL thickness expressed in dimensionless units kcl at volume fractions ip = 1.9 x 10 -2 (circles), ip = 6.5 x 10~ 3 
(squares), ip = 2.9 x 10 -3 (diamonds), and ip = 1.8 x 10~ 3 (triangles), and the correspondent dilute limit extrapolation (stars). 
Left: Transverse motion; right: Longitudinal motion. Curves are drawn as a guide to the eye. 



In Figured we show the sedimentation velocity for a highly charged disk with p = 5, normalized by the sedimentation 
velocity of uncharged disks at the same volume fraction. Although, again, the relevance of the electrokinetic coupling 
in the sedimentation velocity diminishes upon increasing volume fraction, the coupling between electric friction and 
velocity dissipation becomes much more dominant now. The sedimentation velocities decrease by almost 50%, and 
the range of the values of na where appreciable deviations from the uncharged disk behavior is observed is wider than 
in the case of weakly charged disks. Hence, the electrokinetic coupling for disks is much larger than that observed 
for spheres, and is in fact consistent with mobility reductions observed in laponites^l- F° r higher aspect ratios the 
same trends are observed, as shown in Fig. [S] for the particular aspect ratio p — 10. 

In all figures we observe that the reduction in sedimentation velocity for transverse motion is larger than for 
longitudinal motion. For wide double layers the differences may amount to 20%. This effect can be intuitively 
understood in terms of the different forces felt by the electric double layer in the two configurations. For transverse 
sedimentation, most of the diffuse layer is exposed to the flow induced by the sedimenting array of disks. On the 
contrary, for longitudinal motion most of the electric double layer is located in a region where the fluid velocity is 
small and is not subject to large gradients. One would then naively expect that this difference will be enhanced by 
an increase of the surface charge. But, in fact, the relative difference decreases with the charge of the disk. Hence, as 
qualitatively illustrated in Figs. l^land llOl there are non-trivial couplings between the electrostatic restoring force and 
the flow field. From the figures it is clear that for longitudinal sedimentation most of the diffuse layer is in a region of 
smoothly varying velocity, whilst for transverse sedimentation the double-layer is exposed to large velocity gradients. 
More interestingly, by comparing the flow fields past the weakly and the highly charged disk, we observe that the 
effective hydrodynamic shape of the particle becomes more isotropic. Moreover, close to the surface of the particle 
in Fig. EJd, the direction of the flow is reversed with respect to the direction of the bulk flow. This effect can only be 
caused by the different behavior of the electrostatic and the hydrodynamic fields at the edge of the disks, an effect 
we have not analyzed in detail, because much more expensive simulations are required to gain a more quantitative 
understanding of the flow patterns in Fig. 03. Such simulations fall out of the scope of the present paper. 
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FIG. 7: Normalized sedimentation velocity of a charged disk with aspect ratio p = 5, radius a — 5, and valency Z — 100 as 
function of the EDL thickness expressed in dimensionless units Ka at various volume fractions: tp = 3.8 x 10~ 2 (circles), 4.8 x 10~ 3 
(squares), 7.3 x 10~ 4 (pluses), 4.6 x 10~ 4 (stars), and 3.1 x 10~ 4 (diamonds), (a) Transverse sedimentation; (b) longitudinal 
sedimentation. Curves are drawn as a guide to the eye. 
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FIG. 8: Volume-fraction-dependent reduction in the sedimentation velocity of a charged disk with aspect ratio p = 10, radius 
a = 10, and valency Z = 100. The plot on the left refers to a disk sedimenting along its edge, while the right part along its 
symmetry axis. The different curves correspond to volume fractions (p = 1.9 x 10~ 2 (circles), 6.5 x 10~ 3 (squares), 2.9 x 10~ 3 
(diamonds), 1.8 x 10 -3 (triangles), and the corresponding dilute limit extrapolation (dashed line with star symbols). Curves 
are drawn as a guide to the eye. 



VII. SEDIMENTATION VELOCITY OF CHARGED DISKS: SHAPE EFFECTS 



It is not straightforward to quantify the effect of shape variations on the sedimentation of (charged) disks because 
one cannot change the shape of without modifying either the surface charge density or the overall particle charge. 
Then, because the electrostatic field next to a particle is proportional to the surface charge a, a change in the surface 
area will change the electric field surrounding the particle, making it impossible to isolate the effect of shape change. 
On the other hand, keeping a constant by varying the overall particle charge is not a solution either since the reduction 
in sedimentation velocity does depend also on Z [see Eqn. (JSJ]. As a result, we will have to modify both valency 
and volume (to keep the surface area constant) to disentangle charge effects from effects arising from shape changes. 
However, even if we take care of this problem, we can only compare each disk with the corresponding sphere, because 
the two disks we study have different areas. 

In order to focus on shape effects as much as possible, we computed the normalized sedimentation velocity 
Ud{Z)/Ud(0) (with Ud(Z) the sedimentation velocity of an isolated particle with valency Z, and Ud{0) the veloc- 
ity of the same object with Z = 0) with the corresponding normalized sedimentation velocity of a sphere with the 
same valency Z and surface area, U s (Z)/U s (0). 

For weakly charged particles, we can make use of Booth's prediction to analyze the results. To this end, rather than 
studying the scaled velocity directly, we have found fruitful to consider [1 — Ud(Z)/Ud(0)]/Z 2 , which is the coefficient 
C2 [see Eqn. JfJJl] in the case of a sphere. This is a direct measure of the electrokinetic reaction induced by the electric 
double layer. Since we have argued (see Section IIV|I that the charge dependence of disks is the same as the one 
observed for spheres in the Debye-Huckel limit, the previous ratio is a quantitative way of assessing the role of shape 
on the sedimentation velocity. 

In Figures irTT a) and lllf b). we show [1 — Ud(Z)/Ud{0)]/Z 2 = for disks with two different aspect ratios and 
with a small charge, Z — 10, both for transverse and longitudinal sedimentation, cf is expressed in units of A% = 
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FIG. 9: Projection of the fluid velocity field near a sedimenting charged disk on the plane parallel to the axis of revolution of 
the disk. The disk sediment along its edge. 
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FIG. 10: Projection of the fluid velocity field near a sedimenting charged disk on the plane parallel to the axis of revolution of 
the disk. The disk sediment along its symmetry axis. 
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k B Tl B / (727rD?7a 2 ), in such a way that for spheres it reduces to f(na) as predicted by Booth. For weakly charged 
disks and thin double layers, the decrease in velocity does not depend strongly on shape. This is consistent with 
Smoluchowski's theory for electrophoresis 0], which predicts that the electrophoretic velocity of particles with the 
same zeta potential (the electrostatic potential at contact) is independent of the particle shape if na — > oo. However, 
the deviation from this Smoluchowski limit appears confined at narrower double layers for longitudinal motion; hence, 
shape affects significantly the sedimentation velocity of suspended particles. Moreover, in the case of asymmetric 
objects, the orientation of the particle also affects the velocity. For both longitudinal and transverse sedimentation, 
the electrokinetic coupling of a disk is always smaller than the decrease for an equivalent sphere. One can clearly see 
that the decrease in velocity for longitudinal motion is smaller than for transverse motion. 

In the high-charge regime we use the same quantity, , to assess the role of shape, although we know that the Booth 
theory fails in this case. In Figures EI a ) and ll2f b1 we show df, again for two aspect ratios. In the thin diffuse-layer 
limit, our data are consistent with Smoluchowski's theory, and we observe again a departure from the results for a 
sphere upon increasing the width of the electric double layer. The maximum effect is observed for electric-double-layer 
widths of the order of the largest linear dimension of the object, and, again, the decrease for longitudinal motion is 
smaller than for transverse motion. This is consistent with the intuitive picture that the hydrodynamic shape of a 
disk becomes more isotropic upon increasing its charge. 

By comparing Figs. I12f a) and I12f b). the reader might conclude that the reduction in sedimentation velocity is 
higher for the disk with a smaller aspect ratio. However, one should bear in mind that the two disks have the same 
valency and therefore very different surface charges cr p= 5 ~ 3.4 x a p= io- To show how much the surface charge affects 
C2, we show C2 for the disk with p = 10 also at Z = 300. Even though the surface charge of this disk is still lower than 
the surface charge for the other disk, the electrokinetic effect is already more pronounced. 

An illuminating way of displaying the relevance of shape for sedimentation is to consider what the effective Stokes 
radius of a sedimenting disk is. In Fig. ^|we show the effective Stokes radius [R e g = F j (6ttt]U)], where F is the 
magnitude of the external force acting on the disk. For small charges, the effective radius depends weakly on the 
width of the double layer, and is larger for the transverse motion, as can be expected. At high charges the behavior 
is qualitatively different since R c g depends on the width of the double layer for Xd/cl < 2. For larger Xd/cl it tends 
to level off. As the diffuse layer broadens, the effective size that characterizes the sphere and the disk in longitudinal 
motion tend to converge, leading to a same effective shape for wide layers. 

The physical origin of this effect is already implicit in Figs. I§1 andllUI These figures show the velocity fields around 
the disk for both orientations. Different flow fields develop around the sedimenting disk for low and high surface 
charge. The flow profiles look more isotropic for high Z, therefore one might expect that for high Z, the friction 
coefficients of a disk approaches that of a sphere with the same Z . 

VIII. DISCUSSION 

In this paper we have presented simulations of the sedimentation of an array of charged disks. We have treated the 
electrolyte at the Poisson-Boltzmann level, while we have incorporated the relevant hydrodynamic couplings between 
the solvent and the dissolved electrolyte. Using the lattice-Boltzmann method we have modelled highly charged 
colloids and arbitrary na values, which greatly expands the parameter range that can be covered. 

Since no exact analytical expressions exist for the sedimentation velocity of isolated neutral disks and finite thickness, 
we have first checked the performance of our method by validating the sedimentation of neutral disks using approximate 
expressions that become exact for infinitesimally thin disks. Such a computation has provided us with values for the 
sedimentation velocities of uncharged disks, which are needed for the subsequent analysis. 

In order to clarify the role of electrohydrodynamic coupling and the relevance of shape, we have performed a 
systematic study to assess the role of shape, volume fraction, charge, and ionic strength on the sedimentation velocity. 
We find that in the linearized Debye-Hiickel regime, the sedimentation velocity has the same functional dependence 
on volume fraction and surface charge as that for spheres, although with different amplitudes. This deviation should 
be accounted for when using diffusivity measurements of disks to infer the effective charge of colloids and this work 
represent, to our knowledge, the first were the sedimentation velocity is computed systematically. So far experimental 
findings could only be compared with the theory for weakly charged spheres |l ll , which can lead to numerical errors 
in the estimates of their effective charges . At fixed na, we have studied the surface charge dependence of the disk 
sedimentation velocity, from which we have observed that in the high-charge regime, the accumulation of charge near 
the disk surface layer decreases the effect of electrokinetic coupling on the sedimentation velocity, and also shows that 
such accumulation becomes more relevant as the disk becomes more anisotropic. 

We have shown that the geometrical anisotropies of neutral disks are reduced by the presence of the electric 
double layer, especially for highly charged disks. In fact, we have seen that when the double layer is exposed to 
larger velocities, the reduction in sedimentation velocity is larger. Hence, this mechanism tends to generate a more 
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FIG. 11: Normalized reduction in the velocity of a sedimenting disk of charge Z = 10 at infinite dilution both for transverse 
(spheres) and longitudinal (squares) motion as a function of its double layer thickness. The dashed curve corresponds to the 
theoretical prediction for a sedimenting sphere of equal charge and surface area, a) p=5, b) p=10. Lines joining the simulation 
points are drawn as guide to the eye. 
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FIG. 12: Normalized reduction in the velocity of a sedimenting highly charged disk at infinite dilution both for transverse 
(spheres) and longitudinal (squares) motion as a function of its double layer thickness. The dashed curve corresponds to the 
theoretical prediction for a sedimenting sphere of equal charge and surface area, a) p=5, b) p=10. Lines joining the simulation 
points are drawn as guide to the eye. Disks and spheres have a valency Z = 100, which correspond to a surface charge a = 0.45 
for the disk with p — 5 and to a — 0.13 for the disk with p = 10. In Sub-figure (b) we also show the same simulation for a 
disk surface charge a = 0.40 for the disk with p — 10 in transverse (+) and longitudinal (x) sedimentation. Lines joining the 
simulation points are drawn as a guide to the eye. 
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FIG. 13: Effective hydrodynamic radius of a disk normalized by the radius 7? of a sphere with equivalent surface as a function of 
the inverse Debye length in units of the disk radius. Disk with p = 5 in transverse (square) and longitudinal (triangle) motion; 
sphere with same surface (circle). Uncharged objects (drawn line and open symbols), Z = 10 (dashed curve) and Z — 100 
(dotted line). The curves are drawn as a guide to the eye. 



symmetric disk response, as can be effectively characterized in terms of an effective disk radius which becomes less 
sensitive to shape details as charge increases. 
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